QCD string in the baryon 



Yu.S.Kalashnikova|5 

Department of Physics and Astronomy, 
University of Manchester 
Manchester Ml 3 9PL, UK 

Institute of Theoretical and Experimental Physics 
117259, Moscow, Russia 

and 

A.V.Nefediev0 

Institute of Theoretical and Experimental Physics 
117259, Moscow, Russia 

Abstract 

The QCD-motivated constituent string model is extended to con- 
sider the baryon. The system of three quarks propagating in the con- 
fining background field is studied in the Wilson loop approach, and 
the effective action is obtained. The resulting Lagrangian at large 
interquark distances corresponds to the Mercedes Benz string config- 
uration. Assuming the quarks to be heavy enough to allow the adi- 
abatic separation of quark and string junction motion and using the 
hyperspherical expansion for the quark subsystem we write out and 
solve the classical equation of motion for the junction. We quantize 
the motion of the junction and demonstrate that the account of these 
modes leads to the effective "swelling" of the baryon in comparison 
with the standard potential picture. We discuss the effects of the finite 
gluonic correlation length which do not affect the excited states but 
appear to be substantial for the baryonic ground state, reducing the 
"swelling" considerably and leaving room to the short range Coulomb 
force in the three quark system. 
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1 Introduction 



A reasonable model for quark confinement is expected to describe on equal 
footing the simplest quark bound states with zero triality, i.e. qq mesons and 
qqq baryons. While in the naive potential pair-wise picture such uniform 
description is easily achieved by introducing colour factors, in the more so- 
phisticated and more ambitious QCD-motivated approaches one should take 
care of the underlying dynamics of gluonic degrees of freedom when attempt- 
ing to construct the effective model for hadronic constituents. At least the 
non-local gauge invariance should be respected if the quark system is placed 
into a nontrivial QCD background 5^. 

In the mesonic sector the object 'il){y)il){x) becomes gauge-invariant if it is 
supplied with the phase factor (or parallel transporter) $, so that the colour 
singlet qq state is given as 

^{x,y)=r{x)^i{y,x)i,i,{y), (1) 

where 

^i{y,x) = {Pe^V^9 / B^dzX- (2) 

Here a and (3 are the fundamental colour indices, = B^t"' is the back- 
ground gluonic field and P stands to order the Gell-Mann matrices t"" along 
the contour of integration. The only way to form the gauge-invariant state 
of three quarks in the fundamental colour representation is to introduce the 
so-called string junction which serves for coupling three quarks into a colour 
singlet. So the baryonic counterpart of the state (|I]) is presented as 

^(3:1X2X3, Xo) = 

where ^a'/j'^' is the antisymmetric Levi-Civita tensor. The representation 
(H) actually means, that, in contrast to the qq system, the qqq one is to be 
treated as the simplest multiquark system containing a new object — a string 
junction with its own dynamics. 

In the nonrelativistic approach confinement is reduced to the linear po- 
tential acting between the constituents, and the usual assumption about the 
motion of the junction is that the system moves in such a way that the sum 
of distances between the junction and the quarks is minimal (the so-called 
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Torricelli point) [1,2]. This assumption, though being quite natural in the 
potential approach, looks somewhat artificial if we share the common belief 
that QCD at large distances is a kind of string theory and the gluonic degrees 
of freedom manifest themselves as the string ones. 

The latter idea was put on to phenomenological grounds in the flux- 
tube model [2] motivated by the strong-coupling expansion of QCD. The 
ordinary mesons in the fiux-tube model are viewed as qq pairs connected by 
the (nonrelativistic) string in its ground state, and the vibrational modes of 
the string correspond to hybrid mesons. The most recent developments in 
the flux-tube model are given in [3]. 

Another approach is based on the stochastic picture of confinement and 
the Vacuum Background Correlators Method [4]. It is assumed that the 
background field configurations ensure the area law for the Wilson loop op- 
erator giving rise to the straight-line string configuration corresponding to 
the ordinary qq mesons [5], and the perturbations above the background are 
responsible for the transverse motion of the string describing hybrid excita- 
tions. The phenomenological implications of such a model for hybrids were 
discussed in [6]. 

Whatever picture of string confinement one adopts, it is clear that in 
the string-type language the Mercedes Benz string configuration (H) for the 
baryon means that the junction is a special and quite distinguishable point 
of the baryonic string. Its motion might be responsible for special types 
of string excitations which are absent in the simple qq string and reveal 
themselves in systems with more complicated colour structure. This idea was 
first suggested in [2], and the motion of the junction as an independent degree 
of freedom was studied in [7]. Here we continue these studies and present 
more detailed quantitative analysis of the spectrum and wave function of a 
baryon containing four constituents rather than three. 

2 Baryonic area law and effective action for the qqq 
system 

Our starting point is the Vacuum Background Correlators Method. To 
obtain the effective action we write the Green function of the qqq system in 
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the Feynman-Schwinger representation [8] (in Euclidean space): 

roo roo roo r 

G{xiX2Xs,Xo;yiy2y3,yo) = dsi ds2 / ds?, / DziDz2Dzs 

Jo Jo Jo J 



X exp(-/C)- <W >B, (4) 

where 



W = e^,,^-Ari)Kir2WAr,)e-'^'^' , (5) 



3 /.s^ 

/C = ^(m.^Si + - / Zi{T)dT) , 
i=i 4 7o 

with boundary conditions Zi{0) = yi, Zi{si) = Xj. The paths Fj in W run 
over the trajectories of the quarks, and brackets < >b mean averaging over 
background field configurations. 

The Wilson loop operator (|^) can be represented in the form containing 
only closed Wilson loops in the fundamental repesentation by means of the 
relation following from the unimodularity condition for the group SU(3): 

1 = |e„,„,„3^^^''^^^'^'^^(ro)$^:(ro)$^3^(ro) , (6) 

where Fq is an orbitrary open path (the same for all three $'s in (^) con- 
necting points Xq and y^. With inserted unity (H) the Wilson loop (|^) takes 
the form 

W = SpWi SpW2 SpWs - Sp{WiW2) SpWs- 

-SpiW^Wi) SpW2 - Sp{W2W3) SpWi+ (7) 

+ Sp{WiW2W3) + Sp{W3W2Wi) , 

where (Wj)^ is the ordered exponent along the contour formed by the path 
Ti of the i-th quark and Fq (see Fig.l). In what follows Fq will obviously be 
interpreted as the path of the junction. 

ince the relation (^) is the identity for the given background, one may 
rewrite W as 

with Zo{t) running along the contour Fq, introducing explicitly the integra- 
tion over the junction trajectories. 
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To average the Wilson loop configuration (|^ over the background we use 
the cluster expansion method [4] generalized in [9] for the case of several 
Wilson loops or for the loops with self-intersections. It was demonstrated in 
[9] that under the assumption of existence of the finite gluonic correlation 
length Tg for the stochastic ensemble of background fields the generalized area 
law asymptotic can be obtained. For the configuration this generalized 
area law means that the average < W >b is given as 

< W >B~ J Dzo e-'^(^i+^2+S3) ^ 

where a is the string tension in the fundamental representation and Si is 
the minimal area bounded by the contours Fj and Fq In accordance with (|^) 
we integrate over the junction trajectories zq treating in the junction as a 
degree of freedom. The area law (|^) is held for the contours with the average 
size much larger than Tg and it is violated when the contours are nearly 
co-planar. In what follows we neglect such special configurations everywhere 
apart from Section 5 where these effects will be discussed separately. 

The standard approach to the junction motion corresponds to taking only 
the classical trajecory zo{t) in the integral or, equivalently, assuming 
I]f=i Si = min. In the nonrelativistic potential model the latter condition 
is reduced to the assumption that the sum of distances between the quarks 
and the junction is minimal. We would like to note here that for the isolated 
single Wilson loop the area law with S = 5*™" arises in the Vacuum Back- 
round Correlators method if one does not take into account the perturbations 
above the background. However, it follows straightforwardly from the cluster 
expansion [9] that for the configurations under discussion each Wilson loop 
in is actually averaged independently with the result Si = S*™". By the 
way, this result coincides with the one obtained in the expansion (where 
only the first term in survives) , though in the cluster expansion one does 
not apply the Nq —>■ oo limit. Still one needs additional arguments, which 
are not seen at present, to impose the additional constraint J2^^i Si = min. 
The absence of such additional reasons is in fact our main motivation to 
consider the string junction as a baryonic degree of freedom. 

To reduce the four-dimensional dynamics to the three-dimensional one 
we follow the procedure suggested in [5] for the qq system. Namely, we use 
the parametrization 

Zi^ = (r, Ti), Zq^ = (r, fo) (10) 
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synchronizing the proper times of all the particles involved and identifying 
this common proper time with the physical one. After introduction the new 
variables 

/i,(r) = ^iio(r) 0<r<T, (11) 
the Green function (|^) may be rewritten as 

G = I DnDf2Df:iDfoDi^,D^i2DiJ.sexp{-A[{^i}]) , (12) 

where the effective action takes the form 

A= dr 
Jo 

with the surfaces parametrized in the Nambu-Goto form, being the co- 
ordinates of the world surfaces (w,j^ = w-^ = -^)- Assuming, as usual 
reasonable approximation, the straight-line ansatz for the minimal surfaces: 

Wi^ = Zi^{T){l - Pi) + ZQ^{T)j3i , (14) 

with and Zq^ given by the equation (p!OD, we arrive at the Lagrangian 

i-|:g+^+f+-p.|*3.\/rT^). (15) 

where ^ 

k = — [pi X ((1 - (3i)ri + Pifo)] , = fj - fo . 
pi 

Since no time derivatives of the fields /ij enter the Lagrangian, one may 
integrate over /ij in the Green function ( [T^ ) to obtain another representation 
of the Lagrangian 

L = Y. (^miV'l + n' + (yp^ dp,^l + lij , (16) 

that can be considered as the generalization (in Euclidean space) of the qq 
relativistic flux-tube model Lagrangian [10]. 

In contrast to the quark velocities r, the string junction velocity ro enters 
the Lagrangian (|15|) or (|T6|) only via angular velocities Zj. If one neglects 



E 



2/i, 



+ 



Pirt 



— [Wjw: 



l\2 



(13) 
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the terms If under the square roots, then the variable fo becomes a non- 
dynamical one, the integration over fo reduces to taking the extremum in 
To, and the standard requirement [1,2] placing the junction at the Torricelli 
point 

3 

^Pi = min (17) 

is recovered after the Wick rotation to Minkowski space. The particular 
form ( |l5l) of the baryonic Lagrangian with neglect of angular velocities was 
extensively studied in [11], and accounts for angular velocity dependence, 
still within the assumption ([T7|) , was performed in [12]. 

The presence of square root terms in the representation (|1^) makes it im- 
possible to express explicitly the velocities in terms of momenta. To linearize 
the problem we introduce another set of auxiliary fields following again the 
suggestion of [5]. The resulting form of the Green function is given by 

^ = / n Dy^ n Dr,, n n DnOfo exp dTL[{fi, u, 7]}]^ , (18) 

where 



1=1 i=l i=l i=l 



~{\2fii 2 2 Jo 2ui 

+ 11' dPiU.iPifo + (1 - A)rl)' + 1 f dl3iv,rjf{n - ro)' + (19) 
2, JQ Z Jo 

+ d(3iUir]i{fi - fo, (3ifo + (1 - A)^i)^ • 

The initial repesentation is readily recovered by taking the extremal val- 
ues in the fields Vi and 77,. 



3 Effective Hamiltonian of the heavy baryon 

To formulate the Hamiltonian approach starting from the Lagrangian (p!9|) 
one should define the canonically conjugated momenta as 

dL ^ dL 

OTi OTo 
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and perform the Legendre transformation to obtain the Hamilton function 

3 

H = Y1 Pi^i + Po^o ~ L ■ (21) 

i=l 

The Lagrangian (|19D contains only terms which are linear and quadratic 
in velocities, so that the relations (|20| ) may be easily inverted to express 
velocities in terms of momenta. The integration over the auxiliary fields 
/ij, z/j and Tji can be done in the Hamiltonian form of the path integral 
representation as well as in the Lagrangian form (|1^). 

It appears, however, that due to the many-body nature of the problem 
the resulting expression for the Hamilton function is very cumbersome and 
obscuring even before taking extrema in the auxiliary fields and even at the 
classical level. This leads to the delicate problem of appropriate operator 
ordering that would inevitably arise during the quantization procedure. In 
what follows we make several simplifying assumptions. 

First, we consider the case of equal quark masses. Moreover, we will 
restrict ourselves to the lowest hyperspherical harmonic [13] in the quark 
subsystem, the appoximation proved to be successful in few-body systems. 
Anticipating the restriction by the symmetric quark configurations, we as- 
sume that the path integral (|r^) is saturated by the symmetric auxiliary 
fields configurations with 

/ii = /i2 = /i3 = /Lt , z/^ = z/2 = 1/3 = z/ , ni = ri2 = ris = 'n . 

Defining the Jacobi coordinates as 



f 1 = 


73^- 


75^"^ + 


73^2 - 


^^^^ 


r2 = 




75^1 ~ 


73^2- 




r3 = 


73^ + 


V 3^1 " 






ro = 


73^ + 









(22) 



and choosing C, from the condition of mutual cancellation of the terms pro- 
portional to [Rr) we obtain the Lagrangian (in Minkowski space): 

3m^ 3u 3 ,^ 1 ,-2 -2, 1 .9 1,^A9 
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+ /' di3v,i(i - p)({iti + {2(2) + (C - ffn^h - (pii)] , 







where 



12 + Jq dl3v 







Ji = fi+ d(3u{l - (3)\ rris = / d(3u(3{(3 - C) , M = /i + / dl3u (24) 



and p = + 'Ci is the grand hyperspherical radius. 

Another simphfying assumption is to consider the system of heavy quarks 
(m ^ v^)' which allows one to set 

fi = 6 = i? = o, C = o 

everywhere in (pSf) , apart from the quark kinetic term. The integration 
over rj and /i is easily performed, yielding fi = m, the centre of mass motion 
is trivially separated, and the Hamilton function for the heavy baryon in the 
centre-of-mass frame takes the form 



^irip.+ fiff^X^, (25) 

2 Jo Jo 2v 

where p, q and Q are the momenta conjugated to the coordinates ,^1, ,^2 and 
r correspondingly. 

The condition m ^ ^/a allows to treat the system in the adiabatic ap- 
proximation separating the motion of the "fast" string junction and the 
"slow" quark subsystems. Following the adiabatic procedure we rewrite the 
Hamilton function (|25|) as 

2 2 

H = ?,m + ^-^ + H,{p), (26) 
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where the string junction energy Hj{p) takes the form: 

Then the string junction energy Hj as a function of p should be considered 
as the effective potential energy in the quark subsystem, and, similarly, if the 
motion of the junction is quantized, the eigenvalues of the Hamiltonian (|27|) , 
being functions of p, should be considered as the effective adiabatic potentials 
responsible for the quark interaction. 

The Hamilton function (^) of the string junction looks relatively simple, 
but the quantization is not very straightforward: care should be taken of the 
operator ordering in (PTj), and usual Weyl ordering prescription [14] is not of 
the great help here. 

Indeed, usually, when the interaction depends only on the coordinates, 
the ordering ambiguity is resolved in the obvious way. The interaction we 
deal with here is essentially string-type, i.e. nonlocal and velocity-dependent. 
The approximations we have made allow us to localize this feature to a large 
extent, but not to get rid of it: the kinetic term in (p7[) causes difficulties 
in the canonical quantization. Its containing the part proportional to [Qr)"^ 
which explicitly mixes canonically conjugated coordinate and momentum is 
not the whole story — being polynomial in both it can easily be ordered by 
Weyl. The real problem is that the quantization is to be carried out only after 
taking the extremum in u and substituting u = v^xt into the Hamiltonian. 
But the field v enters the effective string mass (p^), and taking extremum 
in V will mix the coordinate and momentum in a severely nonlinear way, 
leaving no hope to get the final answer for the Hamilton operator which is of 
any practical use after Weyl ordering. In what follows we develop the Bohr- 
Zommerfeld quantization procedure, which allows to avoid these difficulties, 
though at the price of losing accuracy. 

4 String junction motion 

The strategy adopted in what follows is to find the eigenvalues of the 
Hamiltonian ( |2^ ) with v as an arbitrary function of (3 and p and then min- 
imize each eigenvalue with respect to so that v = v^xt will depend on 
junction quantum numbers. 
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We start with the classical equation of motion for the junction which 
follows from the Hamilton function (^TD in the spherically symmetric case 



L, = 0: 



z z 



2 



where z = r/p, dots stand for the time derivatives, and the shorthand nota- 
tion 



dp 

10 V 



,2 „2 



1 dfi 



is introduced. The equation of motion (ESf) can be integrated with the result: 



where F{tlj,p) is the elliptic integral of the first kind (see Appendix A for 
the details). The expression (^) defines the function z = z{t) with a and tq 
being the constants of integration. 

Since the velocity i vanishes at 2; = ±a (as it may be seen from eq.(y4.2)) 
the motion is finite and periodic, and the adiabatic invariant 

I= — fQdf=- Q{r)dr 

can be calculated with the help of the eqs.(y4.6), {A.7): 



2 I a z 

I = —nisUjp^Jia) , Jia) = / dz\ — . (30) 

TT JO V 1 + 



2 o X X , \a? - z"^ 

1 + ^2 

The Bohr-Zommerfeld quantization rule for the three-dimensional radial 
motion gives 

/ = 2n + ^, n = 0,l,2,... , (31) 



J(a„) = f- (32) 



so that the expression 



defines as the function of p for each radial quantum number n (^Tj). The 
energy of radial motion is given by 

K(p) = \m,u^p'{al + l) + lf d(3u . (33) 
2 2 JO 



11 



The expressions (0) and (|33| ) solve the problem of quasiclassical quan- 
tization defining the eigenenergies as the functions of p for each n. It is 
impossible in general to write out the analytical expression for the function 
Vn{p)- However some limiting cases are tractable. 

First we note that the quantity a has the meaning of the (dimensionless) 
amplitude for the classical motion, so that the case of small a (a <^ 1) 
corresponds to large p {p ^ r) and vice versa. In the case of asymptotically 
large P (p ^ ^) one expects the string modes to die out. It follows from 
the expression {A.2) for the junction velocity that the limit a <^ 1 means the 
limit f ^ 1, i.e. the "freezing" of the junction around the Torricelli point 
{\n\). Indeed, with the asymptotic (A.S.a) one has 

Vn{p ^ oo) = l-a'p' f'^ + l- f' d(3u + 6Vn{p) , (34) 
/ Jo I' 2, Jo 

where 

5K(p)=^(2n + ^) . 
The extremum in u can be easily taken in (p¥) with the result 



Vn{p oo) = apVs + ^(^2n+^^ , (35) 

so that the spectrum ( |33D becomes the spectrum of the harmonic oscillator, 
which in the leading order in p gives the potential-type linear confinement. 

The actual coefficient of the confinement term in (^) differs slightly 
from the one obtained by taking the lowest hyperspherical harmonic of the 
interaction ([T^). This discrepancy (about 15%) stems from setting all z/j and 
rji equal to each other before quantization and hyperspherical expansion, so 
that our results should be viewed as the variational estimation of the path 
integral (pT^). 

The opposite limiting case of small p is less explicit. One can only consider 
Vn as the function of for a„ 3> 1, what corresponds to small p. Taking 
into account the asymptotic (A.S.b) we have from (|3^) and (^): 

K(an - = i (- + i) + 1 1' d(3u , (36) 

or after taking the extremum in u 

Vn{an^oo) ~ J ^" oo . (37) 

y log an 
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The energy of the string junction grows both for large and small values of 
p, and the position of the minimum of the curve is defined from the conditions 



dp 



P 

V 



Poin) 

^ext 







5v 



P 

V 



Poin) 

^ext 







(38) 



which result in the following expression for Uext'- 



■^ext 



\ 



2(n+|)((a(o)^ 



where a'-"-' 



3J(a(o)) 
2.2 is the root of the equation 
dJ{a 



(a' + 1] 



da 



- 2aJ{a) = . 



(39) 



(40) 



Substituting u^xt from the eq.(^) into (^) and ( PB| ) we obtain the position 
of the minimum of the curve Ki(p) and the value of energy in this point as 



Po{n) 



/ 4n + 3 
cTj(a(o)) 



and 



K(po 



\ 2 



n -\ — 



(a(o)) 



J(a(o)) 



(41) 



(42) 



correspondingly. 

The results for the junction spectrum are obtained by the quasiclassical 
procedure and hence are valid for n ^ 1, apart from the case of large p (^) 
where the quasiclassical spectrum of the harmonic oscillator is exact. The 
quasiclassical results are, however, known to be rather accurate far beyond 
the formal limits of application. So we dare to apply the results to the low- 
lying junction excitations, and even to the ground state. Substituting n = 
into the expressions (ETl) and (H^) we find for ground state 



po{n = 0) 



1.03 



Voipo) ^ 4.8v^ . 



(43) 



It is seen from (^3]) that as the dynamics of the junction is defined by 
the string-type interaction, the typical distances po are determined by the 
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confinement scale Allowing for quark motion will result only in the 

corrections to ( ^31) of order of ^ and we arrive in such a way to a rather 
shocking conclusion: the size of a heavy symmetric baryon is practically 
independent of the quark mass and is essentially defined by the string tension. 
This conclusion holds even for the ground state, and the zero oscillations of 
the string are responsible for the effective "swelling" of the baryon. 

To conclude this Section we briefly comment on the case of Lj 7^ 0. One 
might naively expect that, as soon as the kinetic term for orbital motion in 
the Hamiltonian ( p7D has "normal" quasi-nonrelativistic form, then at least 
for large orbital momenta {Lj ^ 1, Lj ^ rir) the radial motion can be 
neglected, and the resulting energy V{p) will exhibit regular behaviour at 
small p. It is clear, nevertheless, that the case of Lj 7^ is to be considered 
along the same pattern, and zero oscillations in radial motion will again 
result in the effective "swelling" of the system. We have not carried out the 
quantization in this case, mainly because for Lj 7^ there are no a priori 
reasons for restricting by the lowest hyperspherical harmonic in the quark 
subsystem. On the contrary, we expect the total wave function of the baryon 
with nonzero total orbital momentum to acquire considerable contributions 
from configurations with nonzero interquark orbital momenta. 



5 Mass spectrum of heavy baryon 



The Schroedinger equation for the quark subsystem is obtained from the 
Hamiltonian ( p6D with replacing the string junction energy Hj{p) by the 
quantized energy of the junction radial motion Vn{p). If only the lowest 
symmetric hyperspherical harmonic is taken into account, the effective one- 
dimensional equation 



3m — 



2m dp^ 



+ V{p) 



(44) 



with the standard boundary conditions 



Mp = 0) = 



14 



and effective interquark potential 



V{p) 



15/4 
2mp2 



+ Vn{p) 



(45) 



defines the mass spectrum of the system (the necessary details of the hyper- 
spherical expansion are listed in the Appendix B). Note that each adiabatic 
potential Vn{p) gives rise to the whole family of excitations in the quark 



subsystem in accordance with the equation (^). 

The effective potential (^) is shown at Fig. 2 for various values of n and 
quark mass m. The junction eigenenergy V^(p) diverges at small p only 
logarithmically, and the centrifugal barrier dominates the region of small p. 
Nevertheless, since the quarks are heavy (m ^ the minimum of the 
potential V{p) coincides with good accuracy with the minimum of the curve 
Vn{p) given by the eq.(|T]). 

To find the eigenenergies we decompose the effective potential (^51) around 



the point po up to the terms ~ (p — po)^ using the extremal value of v given 



by the expression (|39|) . The resulting spectrum is reduced to the spectrum 
of the harmonic oscillator with the eigenenergies 



en, = V{po)+UJg (uq + 



(46) 



where 



m 



Substituting a^°^ ^2.2 from (|40|) one has 



Vipo) 
V"{po) 



SvrV / 3\ (a(o))2 + l , „ / 4 , 



vrcr 



3/2 



J(a(o))(l + (a(o))2)3 2.7(t3/2 



4(a(o))'^ 2(l + |n) 
with the ground state energy (n = = 0): 

3/2 



1 + ln 



(47) 



ground 



3m 



1 + 1.6^ + 0.9^^'^^ 
m \ m 




(48) 
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The estimate (|i8| ) for the ground state baryonic energy should be con- 
fronted with the eigenenergy derived from the standard approach to the 
junction motion, governed by the interaction ([T7|), with the junction frozen 
at the Torricelh point. In accordance with the discussion of the previous Sec- 
tion, the rehable estimate for this case can be obtained with the first term of 
the interaction (|35|), and instead of ( ^H]) we arrive at (see also [11] for more 
accurate estimates) 



-Aground 3ni 



1.5M^ +0\\^ 
\ m \ \ m 



(49) 



Comparing the spectra (^81) and ( ^9|) one observes that they differ para- 



metrically: the excess over 3m in (^91) is of order of as it should be on 



the dimentional reasons for the linear potential, while in (^S|) it is of order of 
a/o", reflecting the string nature of the interaction. 

Moreover, for the potential case (|17|), (^) the position of the potential 
minimum depends on quark mass — for the ground state one has: 

1 3 

Po = 7 — ^ , (50) 

which tends to zero with increasing quark mass, in contrast to the effec- 
tive "swelling" exhibited by the string approach with po defined from the 



expression (41). 

Before proceeding further with the analysis of the ground state we in- 
troduce short range effects, the first of which is to account for perturba- 
tions above the confining background field. The effect of perturbative gluons 
is twofold. The transverse propagating gluonic excitations give rise to the 
string vibrations and would result in hybrid baryonic excitations similar to 
the mesonic ones [6]. The Coulombic gluons lead to the pair-wise interaction 
of the form 

V-\n,) = \T.CY^^^. (51) 

where the colour factor C is equal to — |, so that the interaction (0) gives 
rise to the attraction between the quarks. 

Averaging the interaction (|5T|) over the lowest hyperspherical harmonic 
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with the help of the eq.(i?.4) we arrive at the effective Coulombic potential 



VSr\p) = -l-. 7 = ^^2.4, (52) 

which should be added to the effective interaction (l45|). 

It appears, however, that for "realistic" values of (in actual calculations 
we take = 0.3) and quark masses the inclusion of the interaction (^) 
leads only to minor numerical changes in the estimates (|48|) for the energy 
and ( 1^ ) for the radius of the ground state, and the higher the state the less 
it is affected. 

There is another short-ranged effect, namely, effect of the finite gluonic 
correlation length Tg. Indeed, lattice calculations give the estimate Tg ^ 
0.2 — 0.3/m [15]. Such a value of Tg is comparable with the typical interquark 
distances, and should be somehow taken into account. 

As was mentioned in Section 2, the area law (H) is violated for the specific 
configurations when all three surfaces or at least two of them are close to 
each other, i.e. the average distance between them is of order or less than Tg. 
The configurations we should be worried about in our spherically symmetric 
baryon are the ones with all the interquark distances less than Tg. 

To handle the problem we first consider the case of coinciding quarks 
trajectories. Assuming the surfaces to be fiat enough we follow the pro- 
cedure suggested in [9] to calculate the average of the complicated Wilson 
loop configurations. Namely, one should find all possible ways to close the 
open contours under consideration, and rewrite the Wilson loop configura- 
tion, using the Fierz identities, as a sum over colour charges in all possible 
irreducible representations of the SU(3) colour group which propagate along 
the coinciding contours, with the coeficients given by the dimensions of the 
representations. Similar results were also derived for the 1+1 gauge theories 
[16] on the basis of completely different formalism of Migdal-Makeenko loop 
equations [17]. 

There is no need, however, to use the whole machinery of [9], because 
the configuration (^ is already one of such representations: the string junc- 
tion is a colour singlet by construction, and configuration (^ in the case 
of coinciding contours corresponds to the colour singlet propagation along 
the common contour. As the string tension is obviously zero in the singlet 
representation (cxo = 0), it means that the string simply vanishes at zero 
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interquark distances: 



< W >B= ^-^^^ = 6 . (53) 

The regimes and (|53|) match smoothly to each other at the interquark 
distances < T^, and we account for this phenomenologically by introducing 
the cutoff function 

f{p) = l-exp (^-^) (54) 

as a multipher at Vn{p) in P5|). Here /(p) is some typical interquark distance, 
and we choose it to be -^p. Such a drastic cutoff of course influences the 
potential curve for small excitation numbers n. Indeed, now the junction 
energy is regular at p = 0. Moreover, there is no minimum for n = 
anymore, while for large n the minimum safely survives the cutoff. Due to 
numerical reasons the minimum appears already for Ucr = 3 at = 0.2 fm. 

The absence of the minimum means that there is no "swelling" for n = 
0, 1 and 2. Such a nonexcited baryon lives in the potential well given to a 
large extent by the Coulombic and hypercentrifugal potentials, and its size 
decreases with the increase of the quark mass. For the standard values of 
quark masses and as this case is qualitatively similar to the case of frozen 
junction [11], as it is shown at Fig.3. Still this apparent similarity is of 
purely numerical origin, and the gap between the energies remains as the 
consequence of zero oscillations. 

To conclude we mention that as the cutoff depends only on the correla- 
tion length, one may find oneself in the amazing situation for large junction 
excitation numbers n > ricr- In the case of large quark masses (numeri- 
cal estimate gives m > rricr ~ ^) with the inclusion of Coulomb force the 
double well potential can appear, and the lowest levels will be the Coulomb 
ones, while another family of levels will populate the (practically undistorted) 
junction curve minimum. 

6 Discussion and outlook 

To summarize, we have demonstrated that allowing the string junction 
to be an independent degree of freedom one does not arrive at the freezing 
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of the junction at the Torricelh point, in contrast to common expectations. 
The eigenmodes of the string junction motion affect the spectrum drastically 
leading to the effective "sweUing" of the heavy baryon in comparison to the 
standard potential model results. The typical size of the heavy baryon does 
not depend on quark masses and is defined by the confinement size The 
standard picture is recovered only if the effects of the finite gluonic correlation 
length are taken into account: in this case the string does not develop itself 
at full scale for the lowest baryonic excitations, and the dynamics is actually 
governed by the Coulomb force in the three-quark system. 

The origin of the effect stems from the special type of baryonic string exci- 
tations, and might be traced to the underlying dynamics of nonperturbative 
QCD. The gluonic degrees of freedom reveal themselves as the junction ones, 
and might be responsible for the spin content of a baryon. 

Unfortunately, (or may be fortunately for the model) the heavy Q-type 
baryons are not observed yet. The excitations wc have described are present 
also in the light quark sector within the string-type confinement model, and 
the effects of "swelling" of the baryon and pushing up the energy will sur- 
vive, though maybe in less pronounced form. For light quarks the adiabatic 
approximation is not valid anymore, and a more adequate approach is to 
consider the four-body system of light constituents with the string junction 
on a par with quarks. Both constituent masses of light quarks and the exci- 
tation energy are of order of ^/a, and one should not expect the quantitative 
features of the four-body system to be dramatically different from the ones 
of the three-body system. As soon as the present quark models define the 
energy of the state up to the additive constant, the increase of energy due 
to junction motion is quite tolerable and might be accounted for by refitting 
the constant. In a similar way, the increase in size which is expected to be 
rather modest for light quarks might be absorbed by re-adjusting the quark 
model parameters. The main effect of inclusion of the junction motion is to 
change the whole pattern of the spectrum due to the presence of the extra 
constituent which carries energy and momentum. 

The most elaborated quark models [18] describe the light baryon spec- 
trum rather well, although some problems remain unsolved, notably the long- 
standing problem of the Roper resonance A^(1440). The mass of this state is 
slightly lower than the expected mass of the first radial excitation, its pho- 
toproduction amplitude [19] does not fit the qqq assignement, and hybrid 
interpretation [20] was invoked to describe this state. It is a challenge for the 
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string-junction model to describe the Roper resonance, together with other 
enfants terribles of the baryonic family like A(1600), as first string-junction 
excitations. 
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Appendix A 

The eg . (^81) can be rewritten in terms of (p = (p' = ^ and z as 



1 



1 + [1 + Z^)^ 

and integrated, giving 

^{z) =uj\l + z^){a^ ~z^) , (A2) 

where a is the constant of integration. It is more convenient to consider z 
varying along the whole real axis from — cxd to +oo, rather then from to 
+00 as it follows from the definition of z, but then while quantizing we are 
to take only the odd levels (see eq.(^Il)). Integrating the eq.{A.2) we arrive 
at the eq.(p9D with 



F(i.,p) = / , "r = /""* , . (A3) 



The period of motion is then calculated using the relation 

F(?/; + 27m,p) = F(7/',p) + 2nK(p) , {AA) 
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where K(p) = F(7r/2,p) is the complete elhptic integral of the first kind, and 
the result is 

T = / K ( , " \ . {A.h) 

To calculate the adiabatic invariant (30) we define the momentum Q{r) 
from the Hamilton function (27) as 

^(^) = -l^'^sT (^.6) 

+ fr 

and, taking into account the relation {A.2), obtain the adiabatic invariant 
(30). The function J [a] entering the eq.(30) can be expressed in terms of 
the complete elhptic integral of the third kind D(p) as 

- vrb° (to) • 

The asymptotic behaviour of the function (A. 7) is given by 

J(a) = -a^ for a < 1 {A.S.a) 

and 

J(a) ^aloga for a > 1 . (A8.6) 



Appendix B 

The hyperspherical expansion [13] implies the decomposition of the three- 
body wave function ^(^1,^2) into a set of spherical harmonics 

6) = E Xk{p)uk{(iIp. 6/p) , {B.l) 

K=0 

where ^^ = ^1+^1) 

(A^- + A^-^)p''uk((i/p, 6/p) = (A^-; + A^-^)Pk = A^Pk = . 



21 



With the decomposition the three-body Schroedinger equation is 

reduced to the infinite set of hyperradial equations for the functions ^k{,p) = 

with 



(5.2) 

Vkk'{p) = J UKV{fi,f2,rz)uK'dVtQ , (S.3) 



where ^(ri,r2,r3) is the potential energy in the three-body system. 

Truncating the set (-B.2) and leaving only the first equation with K 
K' = one arrives at the equation (0). 

The matix element Vqq'^\p) of the pair-wise potential ( ^Tl) 



is calculated from (-B.3) as 

^or(p) = ^^^ — . .. n./L . ; ^ ^ = -7- {BA) 



lo^' d9s tn^9cos^9 J dn^^J dn^jV^^^'in^. f23, rsi] 

Jo^^ dOstn'^ecos^e J dQ^J dQri ' P 



with 7 ^ 
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Figure 2: Effective potentials v= Vn{y^)/ c^^"^ as the functions of x= pa^^"^ 
witli tlie inclusion of hypercentrifugal barrier for m = 1.5GeV (a) and m — 
S.OGey (b). 
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Figure 3: Effective potentials v(x)=K(x)/(T^/^, x= pcr^^'^, with the inclusion 
of hypercentrifugal barrier and Coulomb force for m = 1.5GeV (a) and m = 
b.OGeV (b); upper curve: for the string junction motion with n = and 
Tg = 0.2/m; lower curve: for the frozen junction. 
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